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§ i ABSTRACT 
O 

p ; We present first results of a pilot project aimed at exploiting the potential- 

ity of ground based adaptive optics imaging in the near infrared to determine 
the age of stellar clusters in the Galactic Bulge. We have used a combination 
of high resolution adaptive optics (ESO-VLT NAOS-CONICA) and wide-field 
(ESO-NTT-SOFI) photometry of the metal rich globular cluster NGC 6440 lo- 
cated towards the inner Bulge, to compute a deep color magnitude diagram from 
the tip of the Red Giant Branch down to J ~ 22, two magnitudes below the Main 
Sequence Turn Off (TO). The magnitude difference between the TO level and 
the red Horizontal Branch has been used as an age indicator. It is the first time 
that such a measurement for a bulge globular cluster has been obtained with a 
ground based telescope. From a direct comparison with 47 Tuc and with a set 
■ of theoretical isochrones, we concluded that NGC 6440 is old and likely coeval 

to 47 Tuc. This result adds a new evidence that the Galactic Bulge is ~2 Gyr 
CT) ■ younger at most than the pristine, metal poor population of the Galactic Halo. 
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Subject headings: Globular clusters: individual (NGC 6440); stars: evolution - 



^ ' 1. Introduction 



Galactic globular clusters (GCs) are known since a long time to be the oldest stellar 
population of our Galaxy, hence their accurate age determination is crucial to settle the 
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epoch of the Galaxy formation and a lower limit to the age of the U niverse. A major effort 



has been devoted to obtain reliable ages for the Halo GC s ystem ([Rosenberg et al.l Il999 
Rosenberglliooot ISalaris fc Weis]l2002l ; lDe Angeli et al.ll2005h . 



With the growing awareness that the Bulge GCs are a distinct sub-system, has come 
a new urgency to properly characterize their evolutionary sequences and to determine their 
ages and chemical composition. Moreover, Bulge GCs are key templates of simple stel- 
l ar populations to study the stellar and chemical evolution in the high meta llicity domain 
rtMcWiiliamlll997l : iMatteucci et~al1ll999l : IWvse et al.lbood : iBallero et aDl2007h and for pop- 
ulation synthesis of giant elliptical galaxies. Hence, a proper characterization of the main 
chemical and evolutionary properties of the Galactic Bulge and its comparison with the 
other components of the Galaxy is a fundamental step to un veiling the process of Galaxy 
formation and evolution (IDutra et al.l 120021 ; IWyse et al.l 120001 ) . 



Despite its importance, the Galactic Bulge GC system remained mostly unexplored, 
because of the huge foreground extinction, which severely prevents the usage of optical ob- 
servations. Also, the limited photometric and spectroscopic performances of the previous 
generations of infrared (IR) arrays, prevented any deep, high resolution imaging and spec- 
troscopy to properly survey the Bulge field and GC population. 

Over the last decade a few o ptical and near IR g r ound-based photomet r ic studies of 



19951 : lDavidgell2000[ ). Recently. IValenti et al 



bulge GCs have bee n performed ( Ortolani et al.l Il996l ; iGuarnieri et al.l Il998l ; iFrogel et al. 



(120071 ) presented homogeneous and accurate 



IR color-magnitude diagrams (CMDs) of the evolved sequences of a sample of 24 GCs in the 
Bulge direction. 

However, only a few attempts exist so far to estimate the age of Bulge GCs from the 



onboard HST ( 


Ortolani et al. 


1995 




Heaslev et al. 


2000; 


Ortolani et al. 


2001; 


Zoccali et al. 


2001; 


Cohn et al. 


2002; 


Feltzine & Johnson 


2002 


)■ 



The recent advent of adaptive optics (AO) capabilities on 8m-class telescopes allows un- 
precedent deep ground-based IR photometry down to the TO level in clusters as distant as 
~10 kpc with potential better spatial sampling and coverage than NICMOS onboard HST. 
We performed a pilot project aimed at measuring the age of the Bulge GCs, by taking ad- 
vantage of the ESO-VLT NAOS-CONICA (NACO) imaging facility with AO and wavefront 
sensing in the IR, which is crucial to find natural guide stars bright enough in such a red- 
dened environment. In this Letter, we present the first results of such a project, namely the 
age determination of the GC N GC 6440. This cluster has an iron abundance |Fe/H]= -0.56 
rtOririia. Valenti fc Richll2008h and a global metallicity [M/H]=-0.4 JValenti et alil2007h and 
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it is located in the i nner Bulge at (l,b) = (7.7,3.8) and 8.2 kpc from the Sun. The high redden- 
ing (E(B-V)=1.15, IValenti et al.l 120071 ) along the line of sight yields an extinction of ~3.5 
mag in the visible, and makes this cluster an ideal target to be investigated in the near IR. 



2. Observations and Data Reduction 

The observations have been performed in service mode between August 2005 and Septem- 
ber 2006 as part of the programmes 075.D-0359(A) and 077.D-0316(A) (PI: Origlia), by using 
NACO mounted at the ESO-VLT. We used the S54 camera which offers a relatively large 
55" x 55" FoV with sufficient (0 / .'054pix~ 1 ) spatial resolution to resolve the stars. We selected 
an appropriate field located at «110" South- West from the cluster center, where crowding 
is not too severe but the star density is still large enough to sample a reasonable number of 
stars along the evolved sequences of the CMD (like the Sub Giant Branch). A series of 24 
exposures each one 90 sec-long (DIT=18 sec and NDIT=5 and DIT=10 sec and NDIT=9, 
in the J and H bands, respectively) have been secured. Fig. [TJ shows the NACO J and H 
mosaiced images, with marked the natural guide star used to perform the AO correction. 

We used the final mosaiced images obtained from the reduction p ipeline and computed 



PSF-fitting photometry by using STARFINDER flDiolaiti et ajj|2000j ) specifically developed 
and optimized to obtain accurate photometry in AO-corrected crowded fields. Indeed, 
STARFINDER allows the modeling of the PSF as a function of the distance from the refer- 
ence star with a combination of Gaussian/Moffat functions to simultaneously account for the 
stellar peak and the diffuse seeing halo. More specifically, the complex NACO-PSF in the 
NGC 6440 images has been modeled with a combination of three Gaussians. Two Gaussians 
have been used to properly model the central peak of the PSF (their parameters strongly 
depending on the efficiency of the AO correction), while the third Gaussian parametrizes the 
PSF halo, and depends on the seeing. Since the AO corrections degrade with the distance 
from the reference star, STARFINDER has been optimized to allow the parameters of the 
two Gaussians modeling the central peak to vary as a function of the distance from the ref- 
erence star. All the parameters of the Gaussians have been fine-tuned by properly modeling 
a number of high signal-to-noise ratio, unsaturated stars over the entire science image. Once 
the PSF has been determined, an automatic search of the stars has been performed and the 
PSF fitting procedure executed for all the detected objects. At the end of this procedure, 
a final list of instrumental magnitudes and coordinates for more than 2000 stars has been 
obtained. 



Complementary to the NACO observations of the TO region, we have also acquired 
shallower (i^lmin exposures) J,H,Ks images over a wider (5'x5') field of view at the ESO- 
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NTT using SOFI (programme 081.D-0371(A), PI Valenti) to sample the brightest portion of 
the CMD, namely the Red Giant Branch (RGB) and the Helium clum p. The SOFI J and H 
images have been reduced by using ROMAFOT (jBuonanno et al. 1983 ) and the instrumental 
magnitudes have been calibrated and placed on the absolute astrometric system of 2MASS. 
The brightest, unsaturated stars in the NACO field which are in common with the SOFI 
dataset have thus been used to transform the NACO J and H instrumental magnitudes to 
the SOFI/2MASS photometric system. The average photometric uncertainty turns out to 
be a few hundreths mag along the RGB and ~0.1 mag in the TO region. From a direct 
comparison of our SOFI photometry of the cluster with a control field from 2MASS we 
find that field contamination along the RGB is below 20% in t he NACO field of view. By 
performing a number of simulations with the iRobin et al.l (120031 ) model of stellar population 
synthesis of the Galaxy, we also estimate a similar degree of field contamination in the TO 
region. Hence we conclude that such a degree of field contamination should not appreciably 
affect our estimates of the Horizontal Branch (HB) and TO levels. 



As a template of an old, metal rich GC we took 47 Tuc ([Fe/H]=-0.67) (ICarretta et al. 



2000l ). for which an accurate a bsolute age estimate (11.2 ± 1.1 Gyr) has been recently ob- 



tained by lGratton et al.l (120031 ). As was done in the case of NGC 6440, we u se SOFI shallow 
photometry of the cluster center to sample the RGB (jOriglia et al.l 120071 ) and deep pho- 
tometry of the TO region as obtaine d from archival SO FI J,H images of a field located at 
«2' from the cluster center (see also ISalaris et al.l 120071). The phot ometric analysis of this 
dataset was performed by using ROMAFOT (jBuonanno et al.lll983l ). The sample has been 
astrometrized and photometrically calibrated by using the stars in common with the 2MASS 
catalog. 



3. Results 

Fig. [2] shows the NACO-SOFI J, J - H and H, J - H combined CMDs for NGC 6440. 
Only stars with J < 16.5 in the SOFI sample and within 30" from the reference star in 
the NACO sample are plotted. All the evolutionary sequences from the RGB tip down to 
the MS are sampled. In particular, the CMDs reach two magnitudes below the cluster TO 
and offer one of the cleanest insights of the MS TO region ever published for a cluster in 
the inner Galactic Bulge. These observations show that an AO capability at an 8m class 
telescope is crucial to obtain samples of stars i) ~1 mag fainter and ii) poorly affected by 
incompleteness due to crowding with respect to those obtained from average seeing limited 
observations. These samples allow us to properly measure the TO region even in severely 
obscured star clusters in the Galactic Bulge. 
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The mean ridge line of the RGB/SGB/MS sequence has been determined by applying 
an iterative 3<7-clipping procedure to the observed points and it is overplotted on the CMDs. 

The obtained CM Ds can be used to measure the age of NGC 6440. As has been well 
known for many years (Ubenlll974l ; ISandagdll986l ) the difference in magnitude between the 
HB and the TO is a sensitive function of the cluster age for any stellar population older than 
2 Gyr. In fact in this age range the luminosity of the HB remains nearly constant while 
the MS-TO becomes fainter and fainter with increasing age. The A hb _ MS to parameter, 
defined in different p hotometric bands, has been successfully used in the past to derive 
relative ag es of Halo (jBuonanno et al.l Il989l ; iRosenberg et al.l Il999l ; IChaboyer et al.l Il996l ) 
and Bulge JOrtolani et al.lboOlh GCsT 



In the following we compute this parameter for NGC 6440 and the reference cluster 47 
Tuc, in order to derive the relative ages of these two stellar systems. 

The HB levels in NGC 6440 turn out to be J(HB)=14.70 and H(HB)=13.85, respectively; 
for comparison, 47 Tuc shows J(HB) = 12.5 and H(HB) = 12.05. They have been computed 
by determining the peak of the HB Luminosity Function for stars with J < 15 and (J — 
H) < 0.95. Such a selection allows a safe determination within ±0.05 mag of the mean 
level of the HB red clump by reducing any possible contamination of RGB stars. The MS 
TO points, defined as the bluest point along the mean ridge line of the RGB,SGB,MS are 
located at J(TO) = 18.8 and H(TO)=18.1; for comparison, 47 Tuc shows J(TO)=16.55 and 
H(TO)=16.20. In order to estimate the errors in the TO level determination we follow 



Meissner fc Weissl ( 120061 ). We fitted the ridge line with a fifth-order (J-H)=f(J) polynomial 
and compute the ct(Jto) = cr (f(J))/y(N), where cr(f(J)) is the standard deviation of f(J) 
within a 3er(J — H)to color interval around the TO. The inferred ct(Jto) are 0.07 for NGC 
6440 and 0.06 for 47 Tuc. As an independent estimate we also compute the distribution of 
the Jto values obtained by directly fitting the observed CMD TO region with a fifth-order 
polynomial, by varying the magnitude interval around the TO and the a rejection threshold 
between 2 and 4. The obtained distributions of the TO level determinations are plotted in 
Fig. E The best fit Gaussian a turns out to be 0.08 in NGC 6440 and 0.05 in 47 Tuc. We 
thus conservatively assume a ±0.1 mag error in NGC 6440 and ±0.05 mag in 47 Tuc. On 
the basis of these measurements, the magnitude difference between the HB and the TO in 
the J band of NGC 6440 is AJ(HB - TO) = 4.10 ± 0.11 and the corresponding value in 47 
Tuc is AJ(HB — TO) = 4.05 ± 0.07. These values suggest a similar age for the two clusters. 
The nice similarity between the CMDs of the two clusters is clearly visible in Fig. HI where 
the CMD of 47 Tuc (right panel) has been shifted in color and magnitude according to the 
reddening and distance of NGC 6440, and to match the HB level of the latter (left panel). 
Note that the metallicity of NGC 6440 is slightly higher (^0.1 dex) than that of 47 Tuc. 
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However, inspection of the theoretical isochrones suggests that a 5[Fe/H] w 0.1 translates 
into a 5J(HB — TO) ~ 0.03, only, which is well within the errors. Hence, from this analysis 
we conclude that NGC 6440 is consistent with being coeval to 47 Tuc, that is ~11 Gyr old. 

In order to evaluate the overall uncertainty of the AJ(HB — TO) method in deriving 
absolute ages, we use theoretical isochrones of different ages. Two sets of a-enhanced mod- 
els at [M/H]=-0.65 and [M/H]=-0.35, respectively, to encompass the metallicity of 47 Tuc 
and NGC 6440, and ages ranging from 2 to 13 Gyr, have been retrieved from the BASTI 
database (IPietrinferni et al.l 120041 ) and transformed into the o bservational 2MASS (J,J-H) 
plane following the prescriptions by lOriglia fc Leithererl (120001 ) . The difference between the 
mean value of the HB and the TO level was computed for each isochronal We find that 
in the metallicity range under consideration, a ~0.1mag error in the AJ(HB — TO) value 
translates into a 2 Gyr uncertainty in the 8-11 Gyr age range and >3 Gyr between 10 and 
13 Gyr, since AJ(HB — TO) becomes progressively less and less sensitive to age above 11 
Gyr. Hence we can conclude that NGC-6440 is Gyr, similar to 47 Tuc and other metal 
rich GCs of the Galaxy. 



4. Discussion and Conclusions 



In this letter we have presented the age determination for NGC 6440: it is the first time 
that ground based AO imaging in the near infrared has been used to successfully measure the 
TO region and derive the age of a Bulge GC. The few previ ous studies have been all based 
on space observations. In fact, since the pioneering work by lOrtolani et al.l (119951 ) who first 
measured the TO region of NGC 6553 and NGC 6528 by using HST-WFC2 and suggested a 
near-coeval age for the Galactic bulge and Halo, a number o f subsequent works using H ST and 
the WFPC2 and NICMOS imagers JOrtolani et al.lboQll : iFeltzing fc Johnson! liiioa ) refined 



proper motion decontamination to the CMDs ( Zoccali et al. 


2001; 


Feltzine & Jo 


mson 2002). 


Only one other cluster in the inner Bulge, namely Terzan 5 i 


Orto 


ani et al. 2001 


Cohn et al. 


2002) has been studied with HST-NICMOS. while Heaslev et al. 


(200oh analvzed two other 



GCs observed with HST-WFPC2 in the outer Bulge, namely NGC 6624 an NGC 6637. All 



1 As an independent test of the isochrones, we can also use our me asurement, we can a lso use our 
measurement of the AJ(HB — TO) in 47 Tuc and its robust age estimate by Gratton et al. ( 20031 ) to evaluate 
the reliability of the theoretical AJ(HB — TO) in deriving absolute ages of a stellar system. By using the 
[M/Hl=-0.65 i s ochro ne, we find that 47 Tuc turns out to be «10 Gyr old, which is «1 Gyr younger than in 
Gratton et al.l (|2003l ). In this context, the AJ(HB — TO) value measured in NGC 6440 would give ages sail 



and «10 Gyr, by using the [M/H]=-0.65 and the [M/H]=-0.35 isochrones, respectively. 
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these studies suggest that these clusters are old (>10 Gyr) and likely coeval with 47 Tuc. 



While our and previous studies agree on the fact that Bulge GCs are likely coeval to 
47 Tuc, absolute age estimates range between 11 and 14 Gyr, depending on the adopted 
theoretical isochrones and/or the absolute age for 47 Tuc itself. In this respect, a signifi- 
cant step forward in setting the absolute age of reference clusters has been performed by 
Gratton et al.l (120031 ) who provided accurate distances and absolute ages for three pillar GCs, 
namely NGC 6397, NGC 6752 and 47 Tuc, by using the MS Fitting Method and the Hip- 
parcos distance scale. Their age estimate for 47 Tuc turns out to be 11.2±1 Gyr, w 2.6 Gyr 
younger than the other more met al poor Halo clu s ters. A somewhat younger age (by 1-2 
Gyr) was previously suggested by Rosenberg et al. (jl999 ). who also found slightly younger 
ages for other two metal rich clusters in their sample, namely NGC 6352 and NG C 6838, 
which possibly belong to the thick disk population, like 47 Tuc (see also Figure 2 in iFerraro 
2004h . 



These first set of age determinations seem to suggest a common epoch of formation of 
the tick disk and Bulge and a n overall GC formation process in our Galaxy starting » 13.5 
Gyr ago (IGratton et al.l 120031 ) with the formation of the oldest Halo GCs and ending pa 11 
Gyr ago with the formation of the metal rich objects in the thick disk and Bulge. In this 
scenario, the few younger GCs found in the Halo (like for example Pal 12, Arp 2, Terzan 8 
etc.) could not have formed in situ. In fact, a number of recursive accretions of satellites 
and their GC systems (as e.g. the Sagittarius Dwarf Spheroidal, llbata.Gilmore fc Irwin 
1994T) could have significantly contributed to form the present Halo stellar populations 
Bellazzini. Ferraro fc Ibatal 120031 ; IFerraro! 120041 ) . 



However, a larger sample of both Bulge and thick disk GCs urgently need to be observed 
with the purpose of determining their relative and absolute ages, before drawing firm con- 
clusions. This is something that is well within the capabilities of the current generation of 
instruments, such as the refurbished HST and the improving performances of ground based 
AO imagers. 
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Fig. 1.— NACO J and H band images of the observed field (^50"x50") located «110" 
South- West from the NGC 6440 center. North is up, East is left. The circled star is the 
natural guide star used for AO correction. 
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Fig. 2.— Combined NACO-SOFI (J,J-H) and (H,J-H) CMDs for NGC 6440. Only stars 
with J < 16.5 from the SOFI sample and with r < 30" from the reference star in the NACO 
sample are plotted. The RGB-SGB-MS mean ridge lines are overplotted to the data. The 
HB and TO levels are also marked for reference. 




Fig. 3. — Distribution of the TO level determinations as obtained from fh 400 fitting models 
of the TO region for NGC 6440 (left panel) and 47 Tuc (right panel). The Gaussian best fit 
and the corresponding a value for each distribution is also reported. 
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Fig. 4.— The (J, J - H) CMDs for NGC 6440 (left panel) and 47 Tuc (right panel), as 
shifted in color and magnitude, accordingly to the reddening and distance modulus of NGC 
6440. The horizontal dashed lines mark the HB and TO level of NGC 6440. 



